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Abstract 


We consider the long-time properties of the an obstruction in the Riemann-Hilbert 
approach to one dimensional focusing Nonlinear Schrodinger equation in the semiclas- 
sical limit for a one parameter family of initial conditions. For certain values of the 
parameter a large number of solitons in the system interfere with the ^-function mecha¬ 
nism in the steepest descent to oscillatory Riemann-Hilbert problems. The obstruction 
prevents the Riemann-Hilbert analysis in a region in (x, t) plane. We obtain the long 
time asymptotics of the boundary of the region (obstruction curve). As t —)■ oo the 
obstruction curve has a vertical asymptotes x = ±ln2. The asymptotic analysis is 
supported with numerical results. 
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1 Introduction 

Consider the one dimensional focusing Nonlinear Schrodinger equation (NLS) in the 
semiclassical limit e —)■ 0 


iegt + + 2 = 0, t>0 (1) 

subject to a one parameter family of initial conditions 

q{x,0,e) = (2) 

where real valued A{x) decays fast as x —>■ ±oo and real valued S{x) converges to S± 
as X —)■ ± 00 . 

This is a well known type of problems of finding the leading asymptotic behavior 
of the solution q(x, t, e) in a singular limit. In the case of the semiclassical focusing 
NLS (as opposed to the defocusing semiclassical NLS), the problem is known for mod- 
ulational instability when a smooth initial profile breaks into a seemingly disordered 
structure. 

The first progress in analysis of semiclassical NLS was made by Miller and Kamvissis 
[22j when in numerical studies they observed some order. This has lead to a number 
of results 0 m ED EH [2U] based on the Riemann-Hilbert approach to this completely 
integrable equation. The Riemann-Hilbert approach is based on replacing the nonlinear 
PDE with a pair of linear operators (Lax pair) first introduced by Lax for KdV equation 
[20] and later applied to NLS by Zakharov and Shabat [32]. This reformulates the 
problem for a nonlinear PDE as a scattering/inverse scattering problem for a linear 
operator. So the asymptotic analysis of the NLS becomes an asymptotic analysis of 
the spectral data of a linear operator where the initial data for NLS plays the role of a 
potential. Then the problem is usually further reformulated as a jump (factorization) 
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problem on a contour related to the spectrum of a linear operator - called (oscillatory) 
Riemann-Hilbert problem (RHP). 

Riemann-Hilbert problems are a natural object for the inverse scattering as was 
noted by Shabat [25] who expressed the hardest step - the inverse scattering as a mul¬ 
tiplicative matrix Riemann-Hilbert problem. A (local) RHP we define as the following: 
find a matrix valued function m{z), which is analytic everywhere in the complex plane 
except on an oriented contour E, where the function has a prescribed multiplicative 
matrix jump. Additionally, the function must satisfy a normalization condition at in¬ 
finity. More precise description of the RH approach can be found in nnnnj. A simple 
example of a RHP is the jump matrix to be the identity matrix (’’free” case, no jump). 

To extract the leading contribution, a ^'-function mechanism was introduced by 
Deift, Venakides, Zhou m, applicable to highly oscillatory RHPs. The method is 
based on factoring out contributions until of the remainder RHP has approximately 
constant (upto L 2 correction) jump matrix. Next the ’’model” RHP with the constant 
jumps on finitely many intervals (finite genus) is solved explicitly in terms of Riemann 
theta functions. Then one needs to show that the remainder ’’error” RHP has a small 
(L 2 ) solution. This method can be thought as a nonlinear steepest descent method. 

The Riemann-Hilbert approach to asymptotic analysis has a wide range of applica¬ 
tions to a diverse array of problems including integrable systems (sine-Gordon, Toda 
lattice, (m)KdV, (m)NLS, Benjamin-Ono), combinatorics (longest increasing subse¬ 
quences), Random matrices (GUE, GOE, beta ensembles), and orthogonal polynomials 
(OPRL, discrete polynomials), to name some. 

Eor NTS a number of initial conditions in the semiclassical limit were analyzed 
ig EH ESI [30]. The leading order solution of NTS was found in terms of Riemann theta 
functions with underlying Riemann surfaces of finite genus. Other existing results 
include long time analysis (e = 1) along straight lines | = const in the (x, t) plane 
for several initial data [TIES]. The analysis was similar and lead to a finite genus (0, 
1, and 2) regions in the {x,t) plane. Key ingredients in all these cases were Lax pair 
operators, Riemann-Hilbert problems and the ^-function mechanism. 

Consider the one parameter (^u > 0) family of initial conditions 

g(x,0,e) = -sech (x) gf( 3 ^ 

Even the simplest case = 0 carries many features and difficulties in the analysis. 
A common approach is to approximate the initial data without disturbing the leading 
order of the solution. Moreover, choosing a special sequence e„, —)• 0 leads to purely 
multi-soliton solution which is much simpler for numerical studies. This case was 
analyzed by Lyng, Miller |21j . 

Eor /U > 0, the family of initial conditions Q combines both solitonless initial data 
for /i > 2 as well as radiation in the presence of solitons for 0 < /x < 2. This makes it 
interesting from the point of view of influence of a large number of solitons. 

The solitonless case ^ > 2 has been analyzed completely for all x and t > 0 
values by Tovbis, Venakides, Zhou in |28j . In the semiclassical limit, the leading order 
asymptotics is written in terms of Riemann theta functions with parameters {aj} 
which arise as branchpoints and leading to a Riemann surface. They proved that 
there is a curve t = to{x) (called the first break) in the (x,t) plane such that for 
0 < t < tQ{x) the leading order of the solution q{x,t) depends only on ao{x,t) and 
ai(x, t) = ao{x, t). This can be seen as a genus 0 Riemann surface and the asymptotic 
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Figure 1: Main (bold) and complementary (thin) arcs in the gf-function problem. The exact 
positions of the complementary arcs are not hxed. 


solution has a WKB type approximation. For t > to{x) the leading order depends on 
ao{x,t), a 2 {x,t), ai{x,t) and their complex conjugates (the genus is 2). In the case 
of radiation with solitons (0 < p < 2 ) the previous results were partial: only for finite 
interval of t values (not global in time) |27]. In these studies, there was no information 
on the region/boundary of rigorous applicability. 

In this paper we consider the case 0 < p < 2 with a similar semiclassical approxi¬ 
mation of Q, as it was done in |28j . In these studies the RH approach was completed 
for finite values of t and was not extendable globally for all t > 0. The main obstacle 
for 0 < p < 2 came from a large number of solitons (order of O (^))- These solitons 
correspond to isolated poles of a reflection coefficient of the underlying Lax operator. 
In the semiclassical limit the isolated singularities accumulate and densely fill an in¬ 
terval in the complex (spectral parameter or energy) plane. This adds significantly 
to the difficulty of the asymptotic analysis, which breaks as a leading contributing 
contour coming from analyzing the oscillatory terms in the RHP, collides with these 
accumulated poles and the error estimates become invalid. 

This paper studies the boundary of the region of rigorous applicability of the avail¬ 
able asymptotic result. The boundary t = ts{x) (we call it a singular obstruction curve) 
is a curve in [x, t) plane. We prove that the singular obstruction curve has a vertical 
asymptotes x = Tin2. We find that the rate at which the singular obstruction curve 
approaches these asymptotes isx — ±ln2 = 0(l^) ast—)• -|-oo. We also provide the 
long-time asymptotics of all important quantities including oq, 0 : 2 , 04 . It is conjec¬ 
tured that for |x| > ln 2 the solution maintains genus 2 asymptotics for all t > tQ{x) 
beyond the first break. 

The paper is structured as follows: in section 2 we introduce the main object of 
study - a scalar RHP on ^-function. In section 3 we obtain the long time limit of 
the singular obstruction curve. Section 4 provides numerical evidence supporting the 
asymptotic analysis. In section 5 we discuss the results. Appendix is used for all 
technical asymptotic computations. 
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2 g-function problem 


The ( 7 -function mechanism was introduce(i in HH as a method of extracting the leading 
order by factoring out an unknown function g{z) and setting up conditions to guarantee 
that this function gives the leading order. Usually g{z) is defined through conditions 
on contours in the complex plane. For the semiclassical NLS equation the contour is 
assumed to be a union of so called ’’main” and ’’complementary” arcs on which g[z) 
is defined as the following in the case of genus 2 which is needed for purposes of this 
paper. General setup of any finite genus is similar m- 

1. Main arcs ( 70 , 7 m): 


' g^+g_-f = W 

+ 5 - - / = 0 

9(2(7--/)<0 
^ 9 ( 2 ( 7 + - /) < 0 

2. Complementary arcs ( 7 c): 

f g+- 9- 

\ 9(25 - /) > 0 


on 7m 
on 70 

right from the main arcs 70 , 7 m 
left from the main arcs 70 , 7 m 


on 7 c 

on at least one side from 7 c 


(4) 


(5) 


where U, VF G M. The main arcs 70 , 7m form a branch cut structure which defines g{z) 


In general, the number of main and complementary arcs is determined for each pair 
(x, t), where x and t enter in RHP for 5 -function as parameters through / = f{z, x, t). 
The function / is assumed to be known and it comes from the initial condition 0 
through the logarithm of the reflection coefficient / = ^ log(r). 

In this paper we use / obtained by the semiclassical approximation of the initial 
condition Q as in [28] 


fiz,x,t) = (^1 - 



z+T z—T 

+ In {z + T) + In {z 


T) 


— T tanh ^ ^ — xz — 2tz‘^ ^ ln2, when 9z > 0, 


( 6 ) 


f{z,x,t) = f{z,x,t) when 9z < 0, 

where the branch cuts in the logarithms are chosen as the following: from ^ along the 
real axis to -|-oo, from T to 0 and along the real axis to -l-oo, from —T to 0 and along 
the real axis to — 00 . 


/'(z, X, t) = —^ — In ^ In [z^ — T^) — x — when 9^: > 0 . 


(7) 


In the limit to the real axis from the upper half plane f is 


9/^(2; -t- fO) 


f, z<0 

-|, 0 < z < zGM 

- z> ^ 

2 ’ ^ ^ 2 


so f has a jump on the real axis from Schwarz symmetry. 


( 8 ) 
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Figure 2: Contours of integration in the upper half plane for g{z), h{z) and B{z) in genus 
2 . In the lower half plane the contours are symmetric. Large loop contour 7 and the small 
loop contours 7^, % are shown. For h{z) the point 2; is inside of 7 and outside of 7^, 7c- 
For g{z) the point is outside of 7 and outside of 7^, 7c- For B{z) the point 2; is inside of 
7 and inside of 7^, 7c- As shown, the contours and 2; suit to compute h{z). 


The main and complementary arcs are described by the end points {aj}^^Q: 70 = 
[ai, I] U [f , oo], 7m = [02, 04] U [05, 03], and 7c = [oq, 02] U [03, ai]. Because of the 
Schwartz symmetry of the problem, oi = q-q, 03 = 02, and 05 = 04. Introduce closed 
loops 7o, 7m) and 7c around 70, 7m, and 7c respectively. The orientation of these loop 
contours is clockwise (Fig. [^. The loop % cannot be deformed away from ^ since f{z) 
is not analytic at 2; = ^. We also introduce a clockwise oriented closed loop 7 enclosing 
all the main and complementary arcs together. It is also passing through z = ^. 

The ^r-function have the following expression | 28 ] (Eq.( 3 . 17 ), ( 3 . 18 )) 


9{z) = 


2'Ki 


/(O 


If 




w 


'ym 






( 9 ) 


where point 2 lies outside of a large loop 7 and outside of small loops 70, % and jm- 
The factor R is 


R(^) = \/(^ - ao) (C - ao) (C - “2) - 02) {C - “4) (C - « 4 ), 


where R{^) ~ as ^ +00 and it has branch cuts along the main arcs 70, 7m- The 

genus of the Riemann surface of R{^) is 2 since 70 and 7m form three branch cuts. 
The constants IF and ff are solutions to the system: 


which comes from the requirement g{z) to be analytic at infinity | 28 | . 
The branch points {oq, 0:2, 04} are computed by solving the system 


where 


B{z) 


B{ao) = 0 
B(a2) = 0 
R(a4) = 0 


/(O 






IF 


'ym 


(? - z)R{^) 


d^, 


( 10 ) 


( 11 ) 


( 12 ) 
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with z being inside of the contours 7 , 7 ™, and 7 c (see Fig. [^. 

Additionally the 3 complex branch points satisfy a set of 4 real moment conditions 




d^ = 0, i = 0,l,2,3. 


(13) 


R{0 

These conditions come from expanding g'{z) at infinity into a power series. They are 
necessary but not sufficient to define {ao, a 2 , a^}. The system 0 and formula 
imply (13). 


Introduce a convenient notation 


then 


h{z) = 


Rjz) 

27ri 


fiO 


Hz) = 2g{z) - f{z) 


n 


(14) 


+ 4 (J - 2)R(0 ^ T.,„ ({ - 2)R(0 


w 




where point z lies inside of a large loop 7 and outside of small loops 7 c and 7 ^. 
Then the ^-function conditions ( 0 -®) are rewritten in terms of h{z) in genus 2 : 


1. Main arcs ( 70 , 7 m): 


= 0 

h+ + h_ = 2W 
3/i < 0 
9/r < 0 


on 70 
on 7 m 

right from 70 , 7 m 
left from 70 , 7 m 


(15) 


2. Complementary arcs ( 7 c 


/i_i_ — /i_ = 214 
Sr/i > 0 


on 7 c 

left or right from 7 c 


(16) 


This suggests the visualization: land (9/i > 0), sea (Qh < 0), and sea-shore-lines or 
bridges {Qh = 0). A bridge has < 0 on both sides while a sea-shore-line has Sr/i > 0 
on one side Sr/i < 0 on the other. A main arc can be viewed as a bridge connecting two 
land regions with sea on both sides. A complementary arc is a land path with exact 
position being unimportant as long as Qh > 0 . 

The so called singular obstruction in the procedure occurs when the above assump¬ 
tions that there is contour connecting ^ and —00 consisting of main and complemen¬ 
tary arcs exists is invalid. In particular as we show numerical results in Section 4, this 
scenario occurs at finite values of t for small values of x (|x| < ln 2 ). One of the comple¬ 
mentary arcs collides with the logarithmic branch cut [0,T]. More precise assumption 
on the contours, that there is a path connecting ^ and —00 along which Qh{z) > 0 . 
This condition could fail as shown on Figure]^ for small x values as time increases. We 
call such curve x = Xc{t) in the (x, t) plane the singular obstruction curve defined by 
the condition 

3'/i(T,Xc(t),t) = 0, (17) 

where h{z = T, x, t) is understood in the sense of limit: 


9/i(r, X, t) = O' 


2g{T,x,t) - \imf{z,x,t) 

z^T 


(18) 
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Figure 3 : Schematic illustration of the singular obstruction mechanism in the soli- 

tons+radiation case (0 < /r < 2 ) in the upper half plane. The zero level curve from 0^2 
going up touches the top of the branch cut [ 0 ,T] and closes the passage to connect ao and 
02 (dashed). 

Equation © is an implicit condition for the singular obstruction curve x = Xc{t) 
which we solve asymptotically below. Numerical investigation suggests that this curve 
has a vertical asymptote in the (x, t) plane. In the next section we perform asymptotic 
analysis of the long-time asymptotics of the singular obstruction curve. 


3 Long time asymptotic analysis 


This section contains the core analysis of the paper. It is devoted to solving of the 
equation © in the long time limit. 

Claim 3.1. In genus 2 in the long time limit the branch points ao, 02 and 04 con¬ 
verge to p/ 2 , 0, —p /2 respectively. Convergence of uq p/2 and 0:4 —)■ —p /2 is 
exponentially fast in t as t —)• 00 . 

For p > 2 , the result is proved in [29]. In the case 0 < p < 2, an additional logarith¬ 
mic branch cut [0,r] in the upper half plane appears in equations from which ao, a 2 , 
a 4 are determined. The additional branch cut can be viewed as a small perturbation 
in the limit t ^ 00 and does not affect the leading behavior of the branch points in 
the claim. 


Corollary 3.2. Claim [^/f] implies that the integrals 



no 


(e - z)R{0 


df, 



no 


i^-z)R{0 


dO 



no 


a - z)R{0 


dO 


are exponentially small in time as t —)• 00 . 


no 

n i^-oRio 


(19) 

( 20 ) 


Proof of Corollary 3.2 is exactly the same as in 


Theorem 3.3. Let 0 < p < 2. Assuming the genus 2 region in the {x,t) plane allows 
to send t to infinity, the long time larger solution of the equation 


^h{T, X, t) = 0, 

X > 0 














has the following long time asymptotics 


x{t) = ln2 — 


(8^02 + 11162)62 ^,2 


2 |rp 


— C*3&2 ^ (Q' 262 )) t — y 00 , 


where 0:2 = 02 + ib 2 , T = Cs = ^ 


and 


_ Ai Int 

t 


+ t + i^ + o(ji) 


a2{x,t) = 

62(x, t ) = ^ + ^ + O j 


00 


( 21 ) 


( 22 ) 


with 


Ai = h 


Bi = ^ 


^2 = |ln^ + ^, i?2 = T^ 


, ^3 — 32^, 

or explicitly in terms oft 


where 


and 


2+lnlllli 

R _ ^ I In 2-3; 

-°3 


/ Ini \ 


x(i) = ln 2 + C 2 ^+C 3 ^+ 0 , t-^oo, 


C2 = 


( 8^1 - i?? _ 11 


2 |rp 


C3 = 


8^42 — In Bi 


2 |r|^ 


i-C3)il? = _ 


i 6 |r |2 

// (1 + In 
i 6 |r |2 


Proof: 


(23) 


(24) 


(25) 


(26) 


The proof of the theorem has 4 steps and is presented in the subsections 3.1-3.4. 
The following notation 

J Oj = ?R.{aj), 

I bj = 9(aj) 

will be used throughout the rest of the paper. 


3.1 Simplification of g'{z). 

We start with simplifying the expression for g'{z) for 2 ; values away from 0 and ±/u/2 


9{z) 


m f 

“iTli J ygU-y, 


m 

U-z)R+(0 


dc, 


(27) 


where 2 ; is off the main arcs 70 and 7 m. To overcome the difficulty of R in the denom¬ 
inator taking near-zero values as the cij’s approach the real axis, we transform (27) 
based on the following simple lemma 
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Lemma 3.4. Let C he a closed rectifiable curve in the complex plane. Define 


( 28 ) 


where F{z) is any continuous (on C) function that satisfies the moment condition 


Then 


G{z) = 


ic 


1 


FiOdf = 0. 


z- Zo Jc 

where zo is off the contour C and zq z. 

The proof follows from the simple identity 

1 C- Zo 


(C - zo)F{C) 

i^-z) 


di, 


(29) 


(30) 


f-z {z-zo){f-z) Z-Zo' 


(31) 


which together with the moment condition (|29| ), proves the lemma. 

First we utilize the moment conditions (13) with the contour of integration placed 
along the branchcuts 70 U 7^ 


eno 


d^ = 0, A: = 0,1,2,3. 


These moment conditions can be rewritten in the form 


f M Pmdf = 0 

J70U7m «+(?) ^ 


(32) 


J 7 oU 7 m R+{Q ^ 


where 


f ^^^df = 0 ’ 

J7oU7m rt+(4) ^ 

r miiiidA - 0 

JjoUjm R+(0 U 


MO = “ 2 ), 

^ 2(0 = (^ - «o) if - “ 4 ) , 

A3(0 = if- Oo) if - 0 - 2 ) if - afi) 


(33) 


(34) 


Thus by applying Lemma |3.4| to equation (27), we obtain 

1 Riz) f AsiOfif) 


g'iz) = 


2m Asiz) if - z)R+if) 


df. 


(35) 


Next we deform the contour of integration to the union of oriented arcs in the upper 
half plane. Since / in not analytic on the real axis, oq, a 2 , 04 as points in the complex 
plane, are understood as limits from the upper half plane: oq + iO, 02 + iO, 04 + iO. 
Then 70 U jm is deformed into 

mao) U (ao,ao + ibo) U (02 + ^62,02) U (02,04) U (04,04 + ^64) 
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q;4 


(Xri 


a A ,^ 2 | ao 


20 = f 




as ao 


ai 


Figure 4: Deformed contour Cd- 


and its complex conjugate (with the opposite orientation). We call the new contour 
Cd (see Fig. |^. Then 


9'(2) = 


1 R{z) f A3(0/'(0 


27ri A 3 (z) (^ - z)R{^) 


d^, 


(36) 


and the moment conditions (33) become 


Jcw R(£) U 

fc 


Mono Re - 


Cd R{£) 


d^ = 0 


(37) 


By Claim 3.1, ao and 04 converge to ±/i/2 exponentially fast in the long time 
limit. By Corollary 3.2, the contributions from the intervals (^,ao), (ao,ao + *^o)) 
( 04,04 + ib^) and their conjugates in the contour of integration Cd can be neglected, as 
they are exponentially small as t —)• 00 , and the interval ( 02 , 04 ) can be replaced with 
(“2, — 2 )' Thus (36) reduces to 


9iz) = 


1 R{z) 


Mono 


27ri A^i^z) ^)u(—^,a2)u(a2,02) 2)7?(^) 

which, by recalling f'^ — f'_ = vri on ( 02 , — 2 )’ transforms into 


d^ + exp.small, (38) 


9 [z) = 


R{z) 


I 


^ 3(0 


a2 - z)R{i)^^ ^ -Ki 


1 r Mono 


2A3{z 

Another simplifying observation is 


' 012 


- z)RiO 


d^ 


+ exp.small. (39) 


=^= + exp.small 


MO = 

^(0 -v/(5““2)(€-a2) 

—- _ = + exp.small, 


( 40 ) 


^(0 V(5““2)(€-a2) 

Thus, g'{z) in terms of real and imaginary parts of 02 = 02 + *62 is 


9'{z) = - 


R{z) 

2A^{z) 


{i - M 


'02 


{i-z)\l{i-a2f + hl 


di 
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|•a2—^b2 

+ — 


(e-a2)/'(0 




^iJa2+ib2 02)^61 

Similarly for the last two moment conditions in (|37|) we obtain 

L 


+ exp.small. 


(41) 


1 ra2—ib2 


/'(g) 


, i 1 r 

^“2 V(g-a2)"+fer^ ^ ■Ia2+ife2 y/(^-a2)"+bi 

1 ra2-ib2 (g-ffl2)/'(g) 


+ exp.small = 0 


r-2 _g;i52_Jd , p2-*02 u-«2j/ I _ n 


(42) 


The expression for g'{z) in (41) and for the moment conditions (42) depend only 


on Q !2 = ^2 + ib 2 and not on ao, a^. All the dependence on oq and 0:4 is in the 
exponentially small terms (abbreviated to ”exp.small”). 

3.2 Long time asymptotics of the branch point 02 

To simplify notations we call a = 02 = o(l) and 


We now derive the asymptotics 
6 = 62 = 0 ( 1 ). First we change the variable of integration in the second integral in (42) 
( = by + a 


with 


0^ + o)‘^ + b‘^ 

af + 62 

a 

/(O = 

2 v; 


1 r-i f’jby+a) 


' 771 Jl 




dy + exp. small = 0 


\/?F+I 


1 


(43) 


(44) 


The leading order of /' as / —>• 00 comes from the last term —4t^ since ^ £ [—i, i] 

f'{by + a) =-At{by + a) + 0{l), (45) 

we obtain 


In 


(s+»)"+!-= - (s+») j - ini>+A r +0(1)=0 

^ ^(e+WW-i,+ i/-‘^£^<i9 + O ( 6)=0 


(46) 


By computing the integral [■ * , ^ dy = 0, and expanding the square root for fixed 

* \/y 2 +i 

p > 0 and small a, b. 


In 


2 (f+a) 


+ 0 {!,S))-f/-■^<i!, + O(l )=0 

4tb^ I —i _ 


(f+“)-^r^<i!i+o(;>) = o 


(47) 


Then the integrals are evaluated explicitly f- * ^ ~ —vri, * —^==dy = ^ and 

the logarithm is expanded 


ln5 + 4ta + 0(1) = 0 

f - 2/62 + o(a + 6 ) = 0 ■ 


( 48 ) 
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Solving the second equation for b and plugging the result in the first equation to solve 
for a, we obtain the leading order asymptotics of a 2 = a + ib: 


'>=dt + 0 (l) 


(49) 


As in the pure radiation case, a 2 approaches to 0 on the scale [29]. 

Remark 3.5. The correction term in the second line of (49) is of the order O 
since in the second equation 0{a + b) is in fact 0{a) (see (109), (110) in Appendix). 
The higher order terms are obtained in the appendix: 


a = 5 ?(a 2 ) = + Asi + + O (^) 

b = ^{a2) = ^^,+B2^ + + O (^) . 


This result allows us to compute the long time asymptotics of g'{z). 


3.3 Long time asymptotics of g'{z) 

Let z be such that its distance d{z) from ao and the interval [ 02 , 04] satisfies d{z) » 62 
as t —>■ 00 . First, in (41) we expand in powers of This decomposition uniformly 
holds away from the branch points 02 and 02 , i-e., for \z\ > \a 2 \- Then (41) becomes 


9'{z) = 


1 

vriz 



- 0 ) 


- z)^J{i-af + b'^ 


di 


V(^-a) +62 


+ exp.small. (51) 


Taking into account the last moment condition in (42) in the form 
ca-b^ (^_a)/'(e) 


1 

ni 


-.d^ = - 


^ - a 


J{C-a) +b‘^ 


(e - af + 62 


d^ + exp.small (52) 


we arrive at 


g'iz) = 


(z — o)^ + 62 
2 (z — a) 


(?-a) 




+ A 

^ - Z z ' 


1 

Triz 


ca—bi 


{C-a)f{0 U 




+ E 

k=2 


di 


+ exp. small. 


(53) 


To simplify the expressions, for the rest of this section we reuse (since there is no ao, 

a 4 dependence anymore) notation R{z) = —\J (z — a)^ + 62 and A(z) = z — a in the 
integrals 

" A (0 e 


, 1 R{z) 


R{0 {i-z)z 


di 
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1 R{z) 


ra—bi 


m w 


27riz A{z) Ja+bi 


/'«) ;+E 


k=2 


+ exp. small. 


(54) 


There are two main objects to analyze: and with these ratios 


m)- 


separately by a well known trick 
fn j gn h = {fn - f) J{gn - §) h + {fn - f) j g h + f J{gn - g) h + f j g h. (55) 
Using as guidelines, for t —)■ oo 

R{z) _ af + 62 


A(z) 


z — a 


which we think of as —?• / in (55) and 

m e - a 




1 in 


— 1 as t 


i2 (<.,-f) 


oo, 


— — I as t ^ oo, 


(56) 


(57) 


which plays the role of gn ^ g. 

We do not prove the above statements but rather use them as suggestions in trans¬ 
forming the integrals in (54) 


z)z^^^ 2 


1 fR{z) 


A{z) 


+ 1 


(e - ^)^ 




la U(0 J 




+ 


(? - z)z^^ 2 




1 (R{z) 


1 


A{z 
R{z 


1 R{z) 


z 

f^a—bi 


2'Kiz \A{z) 


+ 1 

+ 1 


_ A(0 

2-Kiz A{z) Jn+bi R{0 


no E 


Kk=2 


U(o Jc-z^ 


Ja+U R{^y Z ^ 

d^ + exp.small. (58) 


For accounting we label these integrals as Ii-e and the terms in the last infinite sum of 
integrals are labeled as k = 2,3,.... We keep terms of the order up to and including 
O (j). See appendix for detailed calculations of the following result: 

' Ii = -5 In (-f - z) + ^ + I In(-z) 0{ay, 

I2 = i2h + 0{aby, 

I^ = ^b^lnb+ (^l^ + yi^y-^ + Oiab), 

U = 0{bHnb) = O , 


h = ^ (2tab‘^ - yfj + 0{ab), Cl = In 
le = 0{taby = O (^) , 

H2 = -^ + o{by, 


2 |T| 


- X, 


(59) 


Hs = 0{tby = O j , 

Hk = 0itb>^+y = 0(^), k = 4,5,.. 
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Recall: b = O and a = O (^) as t —)• oo. 

Now putting together these results 

g'{z) = h+l 2 + h + h + h + h + H 2 + 0{H^) 


( 60 ) 


— /i + O(a^) + + 0{ab‘^) + In 6 + 


l-21n^ ln(l + ^) 


8^2 


+ 


4^2 


+ 0{ab) + 0{h^ In 6 ) H— ^ f 2tab‘^ — ^ + 0{ab) + 0{tab^) — + 0{b^) = (61) 

2:2 Y 4 y 


^ 1 ,2 T 7 2tab^ 3t6^ 

1 /l, , l-2\n^ , ln{l + fi) 


+ 3 I 2'‘ + 


+ 


b^ + 0{ab). 


4 4 

/ 

After substitution of the expression for Ii and Ci the logarithmic terms cancel 

g'{z) = 11 „{ 2 ) - ll„(^+y + ii + i + 2 (a 6 ") 

+ i (IhAlE! + ^+ 4 ) + OM. 


(62) 


(63) 


as t —)■ oo. Note: 0{ab) = O • 


3.4 Long time asymptotics of the singular obstruction 
curve 

In this section we asymptotically solve A/i(T, x,t) =0 which requires asymptotics of 
g{z). By integrating g'{z) to obtain g{z) and using the fact h{z) = 2g{z) — f{z) we 
write ^ 

9{z) = f g'{s)ds +giro) (64) 

Jro 

for some tq G M, ro 7 ^ which implies ’^giro) = 0, that is g{ro) is real. In particular 
we can send tq —)• +00 since g{z) is analytic at infinity. Integrating g' and evaluating 
the limit 

ri 1 . , ,,s ,, 1 

(65) 

( 66 ) 


lim 

ro^cxD 


^ro In ro - ^ (^ro + 0 In (^ro + ^ In tq 


= lim -l(ro + ^)(j^ + 0^^ 


ro^oo 2 


2 ro 


' 4 ’ 


we obtain 


,, 1 , 1 / g\, / g\ g, l f b^ Inb 

g{z) ^-zlnz--[z + -)ln[z + -)+jlnz--[^—+ 2tab- 


1 /I —21n|T| X —ln2\ 2 1 / gb'^ 

'~z V 8 4 j ° ^ V F ”8” 
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( 67 ) 


+ 0{ab) + ^ + g{oo). 


Then for 0 < ^ < 2 and T = i\ 1 — ^ 


%(r) = ^a(T + g 


lnr-ln(r+ I 


+ |r| 


1 /l- 21 n|r| x-ln2» 2 


+ 


1 fh‘^\nb ,2 
b^ + 0 {ah). 


( 68 ) 


and the value of f{T) is computed as the limiting value 


9/(T) = 9mn/(z) = 9('(|-r) 


m 


— + In ^ - T 




+ |r|in(2|r|)-x|r|. (69) 


Next we compute 3‘/i(T) = 2Qg(T) — ^/(T) and after some algebra we arrive at 


3/i(T) = |r|(a;-ln2) + — 


2 / b^ In b 


-\- 2 t(ib^ 


+ 


2 fl-2ln\T\ ^ x-ln 2\^2 

W\ 


Thus the equation 


Qh{T, x,t) = 0 


is asymptotically solved as t —)■ oo 


x{t) = ln2 — 


1 f b'^ In b 


|rp 


+ Atab^ — 


' 4|rp^' ) b^ + Oiab), 


(70) 

( 71 ) 

(72) 


which proves Theorem |3.3| 


Remark 3.6. An interesting observation that in terms of 0(2 = 02 + ^^2 (|50|) 

" Int \ 

,w) 

8 A 2 — In Bi 


,, Int 1 /lnt\ 

a:(t) = In2 + C 2 —+C 3 -+ 0 ( , t^oo, 


where 


C 2 = 


( 8 A 1 - y 


2|rp ’ ” V 2|rp 

To compute C 2 only Ai and i?i are needed and to compute C 3 only A 2 additionally 
required. This may indicate that the terms of the order O (^) and O (^) could be 
combined together in computations. 

Then our conjecture is that the next terms are of the orders O O (^-jj/ 2 '^ 

and possibly O in the singular obstruction curve x{t) long time asymptotics. 

To compute these terms only coefficients: Ai — A 2 , A 3 — A 4 and Bi, B 2 — B^ in the 
asymptotics of 0:2 (t) would be utilized. 


- C3 Bt 


(73) 


(74) 
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X = 0.3, t = 0.6, = 1 



Figure 5: 


Matlab realization of contours of integration for h{z) in genus 2 (see Fig. [^. 




Figure 6 : Contours of integration for h'{z) in genus 2. Large loop contour 7 is on the left. 
Simplified contour of integration 7 ^ for large values of t is on the right. 
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x=0.5, t=1.3, (1=1 X = 0.65, t = 4.0, |i = 1 




Figure 7: Matlab realization of contours of integration 7 for large t values in genus 2 (left). 
Matlab realization of contours of integration 7 ^ for long time computations in genus 2 (right). 


18 
















1 



Figure 8 : '^h(T, x, t, fi) as a function of time t. 


4 Numerical computations 

4.1 Numerical computation the singular obstruction curve 

As discussed in Section 2 the singular obstruction for a fixed x is defined as one of the 
roots of the equation 

A/i(T, X, ts(x)) = 0 (75) 

where the function A/i(T, t) is a nice function of t (see Fig. |^. 

From condition on z being outside of the contour of integration 7 , h{T) is not 
computable directly. Instead 

A/i(T) = A(2<7(r)-/(r)) = 0, (76) 

where f{T) is understood as limit, while function g{z) is analytic at z = T. 

For large t values, the branch point 02 approaches the imaginary axis below the 
point z = T and hits a vertical branch cut [0, T] of function /. Using integration on a 
Riemann surface this event is no special and 02 continues moving on another sheet of 
the Riemann surface without any obstacles. 

4.2 Numerical long time computations 

Fast convergence of ao and 04 to the real axis creates challenges in numerical evaluation 
of gf-function for large values of t. We modify the computations by incorporating our 
asymptotic analysis. 


19 









The branch points uq, a 2 -, and 04 are singularities of the integrands (in the com¬ 
putation of h{z), h'{z), B{z)), and need to be encircled by a contour in the upper 
half plane. This requirement puts the branch points close to the contour of integra¬ 
tion causing falling in accuracy. We utilize here that in the genus 2 region ao and 04 
are exponentially close to respectively as t —>■ 00 and corresponding integrals are 
exponentially small. 

The three complex branch points ao, a 2 , and 04 satisfy the set of 4 moment con¬ 
ditions, which are necessary but not sufficient to compute the a’s: 


^7 ° 

^7 . 


J-y R(0 U 

J 7 


J 7 


d^ = 0 

R{^) ~ 


m) 


(77) 


where 


= Vi^ ~ “0) (? - ao) (C - “2) - 012) - 04) (? - Oi 4 ). (78) 

with the branch cuts chosen along the main arcs: [^,ao], [a 2 ,a 4 ] and their complex 
conjugates. 


By considering linear combinations of the moment conditions (77), the last two of 
these conditions can be written as 


f) = 0 


(79) 


Since ao —^ and a 4 —>■ — we modify the contour of integration accordingly. The 
large loop 7 is reduced to a smaller loop around [—^,a 2 ] and its complex conjugate 
which we call % (see Fig[^. This leads to improvement of speed and stability and 
simplifies the system 


J7r 


J7r 


where 


(?-t)(«+f)/'(«) .7C , 11 n 

^^ ' - 04 -|- exp.small = 0 

+ exp.small = 0 , 

(80) 

(? - 0 “ 2 ) - 02 ) 

(81) 


with the branch cut chosen to connect a 2 and a 2 through — After cancelations, the 


moment conditions (80) look similar to the case of genus 0 with one unknown branch 


point a 2 . All dependence on ao and 04 is in the exponentially small terms. Then we 
approximate a 2 with which satisfies a system of two real equations 




where 


Rlt{ 0 = yj{^- “ 2 '^) - « 2 ^) 


(82) 


( 83 ) 


with the branch cut chosen to connect a^^ and through — 
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One of the key advantages of these long time computations is increased speed in 
exchange of precision in computing 02 • Solving the full system of S-function equations 
0 for ao) ol 2 i and 04 involves computing and inverting a 6x6 matrix of partial 
derivatives for each iteration. While the long time approximations of by solving 
( |82[ ) involve computing and inverting a 2x2 matrix of partial derivatives. 

Using the same contour reduction, we compute U{z) in genus 2 for large t values 


h'lTiz) 


Rlt{z) 

2711 


/'(O 


7r- 


{i - z)RLT{€j 


dC 


(84) 


Long-time computations of '^h{z) in genus 2 as 


QhLT{T) = 3= / h'^rj,{s)ds, (85) 

J ro 

where tq is some real number. 

Long time computations of the singular obstruction curve are based on long time 
computations of a 2 ^(t) and First, {t) is computed which is used to 

evaluate ,t). The long-time approximation of the singular obstruction 

f^^(x) then computed from 

^hiT (T, 02 ^ = 0. (86) 

Using our long time computations of h'^rp[z), and ’^h^Tiz) we avoid the men¬ 
tioned above difficulties with accuracy and improve the speed of the computations in 
the case of and h'j^j,{z). 

The long time approximation ^ of the correct value of 0:2 as a function of time t 
are presented in Figure]^ (left). It shows the time evolution of aoi C( 2 , «4 and for 
X = 1.0 and /r = 1. a 2 and demonstrate similar and converging trajectories. The 
difference convergence tends to zero as t increases. 

The long time approximation tf'^(x) of the singular obstruction curve t = ts{x) for 
/r = 1 is presented in Figure]^ (right). It demonstrates good agreement for t values as 
low as 1 . 


5 Discussion 

5.1 First break 

The first break of the asymptotic solution of NLS Q with the initial conditions Q in 
the semiclassical limit e —>■ 0 was analytically studied by Tovbis, Venakides, and Zhou 
in [28]. They established the mechanism of the first break for the /r > 0 analytically 
as time t evolution process. From t = 0, as t increases the genus changes from 0 to 2 
with a new main arc \a 2 , 04 ] created in the upper half plane. In the present work we 
support their proof numerically. 

Figure |To| illustrates the mechanism of the first break from the point of view of the 
branch points a’s and zero level curves of Sr/i. First, for small t the genus is 0 and there 
is only one branch point in the upper half plane ao. Then, a new pair of branch points 
[ 02 , 04 ] is created, while oq continues to approach ^ under a modified trajectory. 
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Figure 9: Approximation of Q! 2 (t) with cx^'^{t) for x = 1.0, /i = 1 in genus 2 (left). 
Approximation of the singular obstruction ts{x) with [x) for /i = 1 in genus 2 (right). 
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x = 0.4.t= 0.2 ,11 = 1 




Figure 10: Mechanism of the hrst break: t evolution of the zero level curves of for 
X = 0.4, /i = 1. This sequence illustrates the birth of the main arc [a 2 , aj. 
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First break: p, = 0.5 



Figure 11: First breaking curve to{x) (solid) and its small x asymptotics (87) (dashed). 


The asymptotic behavior of the first breaking curve in the large and small x limits 
for p > 2 was established in [28] and 


to{x) = 


+ 


cot - 


2(7+2) ^ ^ ® 
^ _ i In _+ O (-) 


0 


(87) 


oo 


In the soliton+radiation case 0 < p < 2 the above large x asymptotic expression 
produces complex answers. Our conjecture is that the above expression is correct if 
one substitutes T = 0. This conjecture is based on the comparison of our our leading 


order term in the long time limit of a 2 ( 22 ) with the leading term in the expression 
above (1.5) in |29| for p > 2. 

The small x asymptotics in (87) is proved to be valid [55] for 0 < p < 2. Figure 


11 demonstrates agreement of the hrst breaking curve with the small x asymptotic 
formula. 

In this case, the hrst break is a boundary between genus 0 and genus 2 regions. 
From the point of view of genus 2, the hrst break is a singular event of colliding of two 
branch points 02 and 04 with the main arc [a 25 Q: 4 ] reducing to a point. Numerically 


we observed this phenomenon in x, t and ^ evolutions. Figure 13 demonstrates how 
the choice of parameters x and t for the branch points ao, 02 , <+4 evolution correspond 
to the hrst breaking curves. 

The vertical dashed line on Figure [^ (left picture) at x = 0.6 represents the t 


evolution for p = 1 in Fig. 13 (left picture). The line starts above the hrst breaking 


curve for p = 1 so the branch points 02 and 04 do not collide. 


24 
















First break 



First break 



Figure 12: First breaking curves: to{x) for several fi values (left), to(/^) for several x values 
(right). 





Figure 13: Time t (left), space x (middle), and parameter /r (right) evolution of the branch 
points (tto, 0 ( 2 , 0 : 4 ) in genus 2. Collision of 02 and 0:4 corresponds to the hrst break. 
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1 


First break 


First break 




Figure 14: First breaking curves t = to{fi) for several x values. 


The horizontal dashed line on Figure [I^ (left picture) t = 0.8 shows the x evolution 
for // = 1 in Fig. 13 (middle picture). This line intersects the first breaking curve at 


1.5. 


The big black dot at x = 0.8, t = 0.8 represents the parameters for the evolution 
in Figure [T^ (right picture). The dot is located in genus 2 region above all the breaking 
curves which is confirmed by a’s trajectories without collisions for /i values between 
0.05 and 4.0. 

Finally, we look at the first break as a function of the parameter /r. Figure [T^ shows 
no sign of loss of smoothness at /U = 2 which is a critical value for existence of solitons 
in the initial conditions This dependence is investigated in more details in [I]. 


5.2 Singular obstruction 


The mechanism for the singular obstruction is a collision of a branch of zero level 
curves with the logarithmic branch point T (see Fig. 15). This collision closes the 


passage between ao and a 2 around the logarithmic branch cut [0, T] and invalidates 
error estimates. Formally all the expressions h'{z), h{z) are correct as solutions of 


RH problems however the underlying assumptions are not valid: in Figure 15 (right 


picture) there is no path to connect ^ with — oo satisfying Qh > 0 which is necessary 
to guarantee capturing the leading order of the asymptotic solution of NTS. 

It is not clear at this point how to extend the function h{z) beyond the singular 
obstruction. However, the term ’’singular obstruction” symbolizes the difficulties of the 
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x = 0.4, t = 1.583, 


x = 0.4, t = 1,|i=1 x = 0.4,t=1.3, ji = 1 




Figure 15: Mechanism of the singular obstruction: time evolution of the zero level curves 
of for X = 0.4, /i = 1 in genus 2. 


asymptotic analysis rather than drastic changes in the solution. The main difficulty is 
to extend a Riemann surface continuously after a collision of a zero level curve with a 
logarithmic branch point singularity. 

In this paper we confirm numerically existence of the singular obstruction and 
compute two correction terms in the long time limit. 


, In t 1 

X = ln 2 + C 2 — + C 3 - + O 


where 


C2 = 




/ In t 


fi{l + ln 


00 , 


i6|r|2’ i6|r|2 

Comparison of the long time computations of the singular obstruction with the 
asymptotic formula (24) is given in Figure 16 (right picture). 

Based on numerical evidence, our conjecture that the singular obstruction exists 
only for a finite interval [ 0 ,ln 2 ) and the location of the vertical asymptote is indepen¬ 
dent of fj,. 


4|T|^ 


( 88 ) 


(89) 


From the asymptotic analysis of 9/i(T, t) as a function of t in (70), 
^^(T, t) = |r|(x — ln 2 )-|- 0 ( 1 ), t 


00 . 


(90) 


So ^h{T,t) has a horizontal asymptote |r|(x — In 2 ) as t —)• 00 . 

Figure suggests that for x > ln2 function ^h{T,t) has only one root while the 
other root corresponding to the singular obstruction is not present. This conjecture is 
supported by the asymptotics of ^h{T, t) which is asymptotically strictly positive for 
X > In 2 for all t large enough. 

There is no numerical evidence of any other breaks to occur before the singular 
obstruction. It is an open question how to extend the current asymptotics beyond the 
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4.5 


H = 1 



^1 = 1 



Figure 16: Singular obstruction curve t = ts{x) and first breaking curve to{x) for /r = 1 are 
on the left. The long time asymptotics of the singular obstruction curve ( [8^ and the long 
time computations of the singular obstruction curve t = tg'^{x) are on the right. 














singular obstruction curve. More constructively, this is a question of the new genus of 
the Riemann surface after collision of the zero level curve with a logarithmic branch 
point. 

Calculations by Lyng, Miller [2T] in the case n = 0 suggest that after the sec¬ 
ond break the genus is 4. However, they are using a somewhat different approach of 
changing the RHP by adjusting the reflection coefficient. Effectively, they are changing 
approximations of f{z) on the fly: for some region in the (x, t) plane they are using the 
same f{z) as we are, while for other regions they are considering a different approxi¬ 
mation of f{z) (different approximation of the reflection coefficient at later t values) 
and consequently a different RHP. This seems to be equivalent to dealing with two 
/-functions at the same time or dealing with two sheets of the Riemann surface at the 
same time. 

Our analysis is based on approximating f{z) at t = 0 upto order 0{e) and use 
this approximation for all t > 0. R is possible that the leading terms in the used 
approximation of f{z) as e —)■ 0, do not contain information about the second break. 
Thus taking into account correction terms of the order 0(e) and O(e^) is another 
approach to describing the second break. 


6 Appendix 


6.1 Higher order terms of 02 , 62 in Theorem 3.3 


In this section we compute higher order terms in the asymptotics of a 2 = a 2 + ib 2 in 
the long time limit from the couple of moment conditions (42). 


6.1.1 Simplification of the first moment condition (42) 


Consider the first moment condition with the exponentially small terms dropped 


1 


1 


ra—ib 


, -=d/ -| : 


no 


I /-^ 72 

V (^ - o)^ + ^ * 


1 


(91) 


As it was shown above, a = O (^) and b = O j as t —>■ 00 . Our goal is to compute 
the next two terms of asymptotics of a and b. First, we compute the first integral 
explicitly and make a change of variables in the second integral ^ = by + a 


a + ^/(C- a)2 -h b"^^ ^ 




fjby + a) 

VlF+l 


dy 



(92) 


Since both a and b are small, to estimate the integral we decompose ///) in Taylor 
series near ^ = 0 while isolating the term containing t explicitly 


no = -4t^+/'(o) + (/"(o)+4t)e+ +o(e). ( 93 ) 


Next, we incorporate Schwartz reflection symmetry /' (^) = f'(0 and the fact that all 
/(O(o) are real for k > 2, then 

f'{by+a) = -4:t(by+a)+Ci+iC2sign (^^^+C3{by+a)+Ci{by+a)‘^+0({by+a)^), (94) 
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(95) 


where the constants Cj are obtained from (|44|) as limiting values 


So 


= lim + ie)). 

£-S>0+ 


Cl = In — X 


TT 

02 - 2 
Cs = - 

■5 


(96) 


5+“l 




After plugging (|94|) and (|96|) back in (|92|) one obtains 
In 

1 / 

+ 


1 r^-U{by + a) 


+1 

1 /■“* C*! + iC 2 sign (f) + C'i{hy + a) + Ci{hy + af + 0{{by + a)^) 


dy 


dy = 0, (97) 


Ji Vv^ + i 

where after expanding all the terms, using symmetries in the integrals, and only keeping 


terms upto order 0(6^) = O (]-), we 
In 


see 


X^+a) 

which simplifies to 

In 

and its final form is 


+ 0 ( 6 =)) + 1 r + 0 ( 6 =) = 0 , 

m Ji / n . , 


+1 


2(f+ a) 


- (-4to + Ci + Cga) + 0{b^) = 0 


In 6 + 4ta — Ci — ln//+( - C 3 ] a + 0(6^) = 0. 


Finally, we plug in values of Ci and C 3 from (96) 


Inb + Ata + x- ln2|r| - — + 0(62) ^ 


(98) 


(99) 


( 100 ) 


( 101 ) 


In order to extract the asymptotics of a and 6 from this equation we need to conple it 
with the other moment condition. 


6.1.2 Simplification of the second moment condition (42) 

Consider 


/■ 


^ - a 


1 


ra—ib 




! - -d^ H : , - 


di = 0 . 


( 102 ) 


In a similar manner as for the other moment condition, we evaluate the first integral 
and make a change of variables ^ = 6y + a in the second integral 


^ + a) +62-6+A 

2 / 7^^ 


1 


■* byfjby + a) 

\/y2 + 1 


dy = 0. 


(103) 
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Next, we compute the integral by expanding f'{C) into Taylor series (94) 

/■-' byfiby + a 


-dy 


(104) 


Ji 1 

by [-At{by + a) + Ci + iC2sign (f) + 6*3(6^ + a) + Ci{by + g)^ + 0 {{by + af)] 

\/y‘^ + 1 

(105) 

and using the symmetry in the integral which make many of the terms to disappear 


i: 


= {-Atb^ + Csb^) 


+ 1 


dy + 2 C 2 bi 


+ 1 


dy + 0 {ab^ 


= {-Atb^ + C^b^) y + 2C2hi + O {ab^) 


Then the equation ( 103| ) reads 

\/(f+a)T6=-6-2tt'^ + ^ + ^ + O(»l,2)=0 

and after expanding the square root 


/f 

2 




(106) 

(107) 

(108) 

(109) 


This equation is an extended version of the second equation in (48). After substituting 


expressions for C 2 and Cs from (96) we see that the term of order the 0{b) vanishes 

( 110 ) 


— — 2 tb^ + a + + O iab^) = 0. 

2 y 


262 


6.1.3 Asymptotically solving the system of moment conditions (42) 

We want to solve the system of two asymptotic equations for a and 6 


In6 + 4to + x-ln2|r| - ^ + 0(62) ^ q 
^ — 2 tb^ + a + ^ + O (ab^) = 0, 


( 111 ) 


where as we established in (49) the leading order solutions are 

/ «=W + o(i) 

1 '> = \/S + o{i). 

We find the correction terms by writing 

{ _ InJ. £1 Int 

“ ~ 84 t 

6 — . 4 - 

with functions ei, 61 to be determined. Then 


( 112 ) 


(113) 


ln6 = Int + In + 0((5i) 

4,2 _ ii I vT'?i , 4 
^ — 4 t ^ 4 4 ’ 


(114) 


dy. 
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which we substitute into ( 111 ) and obtain 


I — ^ Int + In y f + 0(5i) + 5 Int + 4ei hit + x — ln 2 |T| + 0 {a) = 0 
\ f f - 2^2 + ^ + ^ + O ( 62 ) = 0 . 

This leads to the following asymptotic equations 


(115) 


deilnt + lny/f+x-ln2|r| +0(6i) + 0(o) = 0 
-2^6i + ^ + O + O (62) + 0(62) ^ 0 ^ 


(116) 


which we solve for ei and 6 i 


_ ip^/f+3:-in2|r| 

~ 41nt Ul 

+o(i) 


6 i = 




and thus 


, , ln^lll+ln 2 —X 

“ ~ 8t 4t 


+o{'m 


b= i/* + 


Int 


+ o 


(117) 


(118) 


4t ^ 

In a similar manner we compute the next order terms for a and 6 by introducing 
correction terms with unknown functions 62 and 62 


_Int I A 2 I £2 Int 

“ “ 8t t “T *2 

h — /it 4 - ^2 Inf I (52 
" V 4t t3/2 f3/2 ■ 


(119) 


where 7 I 2 = ^ In ^ 


and B 2 = 16 ^ then system (111) becomes 


2 B 2 Inf I 4£2 Inf _ ln_f p, /l'\ _ p. 
^f f 2fit ^ ^ \ t) ~ ^ 

-^ + t + s + o{¥)=o. 


( 120 ) 


which we solve for 62 and 62 ^ 


^2 - ^ + 0 (i^) 


62 = 


^7^ 


+ o(!y). 


( 121 ) 


This provides the three leading terms of the asymptotics for a and 6 

2|T| 


a = 


4 v7 ^ 16Vif f3/2 
which completes the asymptotics of 02 . 


, , , In-LA-j-fin 2 —x) , 01 + / 1 \ 

llnf, ^ ' 1 , 

8~ ^ 4 f + ^ 

2+ln^ + (ln2-x) ^ 


^= 4/^45 + 


1 Int 


+ 


8V7 


f3/2 


+ 0 


(^) 


( 122 ) 
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6.2 Technical Lemma: asymptotics of an integral 

Lemma 6.1. Let f{s) be twice continuously differentiable on [0,p] with some p > 0 
and f"'{s) exists and bounded in a small neighborhood of s = 0 then 


j f{s) (^V + 62 _ (Ig = 


/(O) , /(0)ln2 f{p)lnp 1 fP , 


+ + 


+ 


f (s) In sds'^b^+ 0{b^), 6^0+. (123) 


(124) 


4 ' 2 ' 2 2 

We just outline the proof; 

1. Split the original interval into two subintervals: small neighborhood near zero 
and the rest 

[0,p] = [0,(5] U [S,p], 

for some ( 5 — 7 - 0 '’“. 

2. Show 

b^ 


(126) 

(127) 



/I In 2 

ln(5^ 

1,0 ^ ( 

V4 + ^ 

H—^ 

2 y 

)b^ + Oi 


J s"^ + 6 ^ — s'j ds = O if?5^) . 


3. Show 


j f{s) s2 + 62 - S^ds = /(O) 


6^ In 6 /1 In 2 In (5 

—^-1" 1 T H—^—I— 


4 2 


+ f'{o)^^ + o{b^)+o{b^s^). 


(128) 


4. Show 


f{s) (^Vs^ + - s) ds = 6 ^ ® ds + O (^ 6 ^ 


5. Show 


^ f{p)^T^P /(0)ln(5 /'(0)(5 


Is 2s 2 

1 fP 
~ 2 


2 2 
r f'{s)lns ds + 0{S^), 

Jo 

m 


Is s" 


-ds = O ( p 1 . 


(130) 

(131) 


6 . Show 


= b 


CP 

/ /(s) s^ + 6^ — s^ ds (132) 

, ^/(p)i„p_ f'M In»*) + + (?())). 


V 2 

7. Set (5 = Vb. 
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6.3 Asymptotics of key integrals 

Here we derive the asymptotics of the integrals in 


11 = -\\n (-f - ^) + 4 i + 5 ln(- 2 ;) + 0{a?), 

1 2 = + 0{ab'^), 

h = 4 ^ 6 ^ 1 n 6 + b"^ + 0{ab), 

U = 0{bHnb) = 0Cf), 

h = ji (2tab‘^ - + 0{ab), Ci = In ^ - x, 

h = 0{tab^) = O , 

H2 = -^ + 0{b^), 

Ffc = 0 (t 6 "+ 2 ) = O , k = 3A,... 

in the long time limit t —)• 00 (a —>■ 0, 6 —)• 0). 


(133) 


6.3.1 Asymptotics of Ji, I 2 

Let \z\ > |a|. Consider 

— 


u 


1 1 

7-^ 

^ — z z 


and then isolate the leading order 0 ( 1 ) 




'2 


Ml 

4z 2 


—di 


(134) 

ZjZ 



2 

a 

(135) 

>] + 5‘“('-i) 

a 

^Yz 

(136) 

ln(—z) + O(a^). 


(137) 


The correction term is of the order 0{o?) = O 
For the second integral 


InO 


h 


^ _1 (RA) 


2 VA(z) 

Consider the front factor first 


+ 1 


{i-z)z \k{z) 


(138) 


for \z — a\ > \b\ 


for \z\ > la 


Rjz) 

A{z) 


+ 1 = 1 - 



62 

{z — a )2 


+ exp.small 


= [1 + 0{a)f + 0 ( 6 "^) 


(139) 

(140) 

(141) 
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— + ^(®)] + 0 {b'^) — + 0 {ab^). 


Thus, since /i = 0 ( 1 ) 


I2 — 2^^^ 0 {ab^). 


(142) 

(143) 


6.3.2 Asymptotics of I 4 

Consider 


(m 




1 

2 


V^(6 / iC-z)z 

C - a 




- a)2 + 62 
Next, we perform integration by parts 


+ 1 I I 7 I ) d^- 

X — z z 


(144) 

(145) 



(146) 


{s + a — 2:)2 


and we keep only terms of the order up to 0 ( 6 ^) = O (7), 0 {ab) = O 

IJ 62 


ds 

(147) 


4 (f+ z)z 


\j{^ + a) +62 + (f + a) 




(s — a + z)2 


4 (^ + 2 ;) 2 ; 2 Jq {s - a + z)^ 


(149) 


In the integral a enters as a part of the upper limit of integration and as a part of the 
integrand. Introduce a notation 




1 f2+'^^MT¥- 


then 


(s — a + 2:)2 
J(a) = J( 0 ) +j;( 0 ) a + 0 {a^) 


ds 


1 [2 a/s^ + 62 — 


ds + a 


2./0 (s + ^)2 I 2(f + z)^ 

1 fa a/s^ + 62 — S 


2 ) “ 2 _ 1^ Vs2 + 62 - S 


(s + z)^ 


ds + 0 {ar) — a 


2 do (s + 2^)^ 

72 a/s 2 + 62 — s 
/o (s + 2;)3 


ds. 


(150) 

(151) 

+ 0(a2) 

(152) 

(153) 
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Both integrals have a similar structure f{s) ds which was studied in 

the previous section. From Corollary |6.1t the first integral has asymptotics 


r/i/2 


2 (s + z) 


+ 


1 — 2 In ^ 1 


+ 


^ (V+ 62 _ s) ds = ^b^lnb 


1 


+ 


8^2 2nz 2/i(z + 0 4z2 

and the second integral in ( 153[ ) has the order 

^ ds = 0 ( 6 ^ In 6 ). 

Jo (s + z)^ 


6 ^+ 0(6^), 6^0 (154) 


Thus, returning to (149) 


+ 


1 — 2 In ^ 1 


l3 = 


+ ■ 


-r-h 0 {ab) + —^b'^lnb 

4(f + z)z 4^2 


1 


+ 




(155) 


6 ^ + 0(6^) + 0(afe^ In 6 ) (156) 


^ 8^2 2 fj,z 2/i (z + ^) 42;2 ^ 

and since 0{ab) = O 0{ab‘^\iib) = O (^) and 0{b^) = O 

62 62 1 


+ 


2^(^+2;) 2^z 42;2 


+ ^ 6 ^ In 6 


+ 


l-21n/r 1 1 ln(l + ^)\ „ ,, 

+ —-U7V H- \ o & + 0(a6). 


8z2 


2/iz 2;U (z + 4^2 


So 


I3 = ^62ln6 + 


' l- 21 n/r ln(l + ^) ' 
8^2 42;2 


62 + 0(a6). 


Consider 


, _ 1 (’iOO , 1 


U® 




Similar to our calculations in Step 1 

= -(-T+0(a®)/3 

and since I 3 = 0(62 In 

/4 = 0 ( 6 ^ In 6 ). 


(157) 

(158) 

(159) 

(160) 

(161) 
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6.3.3 Asymptotics of Js, Jg 

Consider 


h — 


1 


1 




2 'Kiz‘^ 


a+u — a)^ + 6^ 




changing variables ^ = 6y + a 


1 


2 Txiz^ 


byjby + a)f{by + 

+ 1 


dy 


decomposing /(^) into powers of ^ as in (94) 

/■■* + a) 


1 

2 'Kiz'^ 


-4:t{by + a) + Cl + iC 2 sign 


\/y^ + 1 

+C' 3 (^d + n) + C/i{by + a)^ + O {{by + n)^)] dy, 


where (see (96)) 


Cl = In — X 

C TT 
2 - 2 

cs = i, 


leading to 


(162) 

(163) 


(164) 


(165) 


(166) 


1 f ' {b‘^y‘^ + aby) 

2vri2;2 ^y2 _g i 


—Athy — Ata + Ci + iC 2 sign ( v ) + 0 {h) dy 


(167) 


by the symmetry argument all odd powers of y vanish 


^ ^ +l + Cib^y‘^ + 0 {ab)] dy 


J. ^y 2 

{—Stab"^ + Ci6^) 




2Triz‘^ 




dy + 0 {ab). 


Thus since A * T— dy = ^ 


/s — —^ ( 2 tab + 




+ 0(a6). 


Consider 


l 6 = - 


1 /d?(z) 


27riz VA(z) 
R{z 


+ 1 


f^a—bi 


MO ^ 


la+bi MO^ 


A{z) 


+ 1 Is = O{b0O{tab0 = o 


Int 


(168) 


(169) 


(170) 


(171) 

(172) 
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6.3.4 Asymptotics of 

Consider 


Ho = - 


1 R{z) 


2 mz K{z) Ja+bi R{0 




1 pa —01 




_ (e - af + 62 

with the change of variables ^ = by + a 

_ 1 + 0(6^) byjby + affjby + a) 

2vri2;3 ^y 2 i 

similar to computation of Is 

1 + 0(62) r-i (53y3 + 2a62y2 + a?})y) 


dy 


(173) 

(174) 

(175) 


2 'Kiz^ 

1 + 0(62) i-i 

2 'Kiz^ I 


\/l/^ + 1 


[—4t6y — 4ta + 0(1)] dy (176) 


\/y^ + 1 


[-4t6^y^ + 0(ta6=^) - 4ta63y3 + 0(ta262) + 0(6^)] dy. 

(177) 

By the symmetry of the integral the y^ term vanishes as well as the term 0{tab^), 
leaving 


/f2-(l + 0(6=))|h +0(1,3) 


2 tb^ /•-* y^ 


niz^ 


\/y^ + 1 


\/y^ +1 

dy + 0(6^) + 0(f6®) 


(178) 

(179) 


with a table integral f. * ^ 
So 


^dy = 


3t6^ 

42:^ 


Consider for /c = 3,4, • • • 


Hu = 


1 d?(2) 


similarly to H 2 


-I pa —01 


■ + 0(63). 


(180) 

“A(0 /O 

k 

) di 

(181) 

-{(-a) 

1 

=notdi 

(182) 


(^ - a)2 + 62 


with the change of variables ^ = by + a 



_ 1 + 0(62) r * 6y(6y+ a)*'/'(6y + a) 

27ri2*^+i Ji y/y 2 + 1 ^ 

(183) 


1 + 0(62) r-i 0{b)0{9)0{tb) 

2 mz^+^ Ji \/y2 + l ^ 

(184) 

Thus 

Hu = 0{tb^+^) = 0{b’^) = 0{b^). 

(185) 
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6.4 Numerical evaluations 


6.4.1 Numerical evaluation of h'{z) and h{z) 


The main idea we utilize is implementing integration of functions on Riemann surfaces 
rather than in the complex plane. Such approach together with numerical contour 
deformations has allowed to avoid expensive computations of the main arcs as a prepa¬ 
ration for any single computations involving /i-function. We were able to continuously 
track the branch points a 2 in genus 2 beyond colliding with the branch cut [ 0 , T] to 
another sheet of a Riemann surface. This leads to easier long time computations and 
allows to observe the singular obstruction. 

To compute h! [z) we assume that all main and complementary arcs are in series 
configuration, and the positions of the branch points a’s are known (see below), then 


h'{z) 


Rjz) f no 

2vri (? - z)R{0 


where the point z is 


inside of the large loop 7 (Fig. and where 


fiO = -y “ ( 2 “ ^ -T^) - X- At^, when > 0. 


(186) 


(187) 


Under the same assumptions, we evaluate the function h{z) directly rather than 
integrating the derivative h'{z). 


h{z) = 


Rjz) 

27ri 


fiO 


L./7 


- z)R{C) 


d^ + 


n 


7c 


{(-z)R{C) 


-|- 


W 


'7m 


(C - ^)R(0 


di 


where point z lies inside of a large loop 7 and outside of small loops 7 c and 7 ^ (Fig. 

©• 

The constants and W are solutions of the linear system 


^7. 4 )^^ 4 4)^0 V ^ ^ 


(188) 


The code to evaluate h{z) could be viewed as a tool to support and track evolution 
of RH contours. It provides valuable insight into the behavior of the branch points 
and the zero level curves of 71/i. This tool allows one to identify which one of several 
possible scenarios of level curve evolution does occur. 


Remark 6.2. The only principal difference between evaluating h{z) and g{z) is a re¬ 
quirement for the point z to be located either inside contour 7 (for h{z)), or outside 7 
(for g{z)). 

Keeping in mind simple relation between h = 2g — f, ii is easy to switch between 
these two functions. For example, for distant z it is more efficient to use g. 


6.4.2 Numerical evaluation of o’s 

We compute the branch points (a 0 )« 2 ,« 4 ) in genus 2 by solving the system 

B{ao) = 0 
B{a 2 ) = 0 , 

B{a/i) = 0 


(189) 
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where 


B{z) = f 


/(O 


n 




W 


'Ym 


- ^m) 


de 


(190) 


Up to a constant, B{z) is h{z) without the R factor in front of the integrals. 

Remark 6.3. We stress that B{z) is in fact a function of a’s and a’s through the R 
factor in the integrals, that is B{z) = B{z,a). While B{z,a) is an analytic function 
of z, B{a,a) is a non-analytic function of a which depends both on a and a through 


R{^) in the denominators. So we treat the system (189) as a real 6x6 system and solve 
it iteratively. 


6.4.3 Numerical computations of the first breaking curve 

Computations of the first break for fi = 0 was done by Lyng and Miller |21j . We 
computed the first break for /x > 0. The first break was also observed as a singular 
event in x, t and even /x evolution of the zero level curves of 9/i. 

We treat the first breaking curve t = to{x) as a function of x. For fixed /a and x 
we are looking for a pair {zo,t) which satisfies the system of one complex and one real 
equations 

h'{zo,t) = 0 

‘^h{zo,t) = 0 , 


(191) 


where ^ 0 and t > 0. The Jacobian of this system is singular. We start with 
expending h{z) in powers of 2 ; — 


(j 

h{z) = A + B{z - zo) + —{z - zof + 


(192) 


where A = h{zo), B = h'{zo) and C = h"{zo). Then the system is approximated as 

B + C{z — zq) = Q 


[A + B{z - zo) + ^{z - zo)^] = 0, 

Solving the first equation and substituting into the second equation leads to 


A-^ 

2C 


= 0 . 


(193) 


(194) 


Thus in terms of the function h the system (191) is replaced with 


f h'( 

zo,to) = 0 



h{zoRo) - 

(h'{zQ,to))'^ 

2h"{zo,to) 


= 0 


This system is solved iteratively where zq and to ^-re updated in turns 




(1) 


hi) 


.(2) 


(195) 


(196) 


We use the first equation in (195) to update zq and the second equation to update to- 
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6.4.4 Numerical computations of the singular obstruction 

From the numerical point of view, the singular obstruction curve is a solution of a 
scalar equation 

^h{T,x,t) = 0 (197) 

for either t = tc{x) or x = Xc{t). 
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